Clinical Implementational and Site-Specific Workflows for a 1.5T MR-Linac.

MR-guided adaptive radiotherapy (MRgART) provides opportunities to benefit patients through enhanced use of advanced imaging during treatment for many patients with various cancer treatment sites. This novel technology presents many new challenges which vary based on anatomic treatment location, technique, and potential changes of both tumor and normal tissue during treatment. When introducing new treatment sites, considerations regarding appropriate patient selection, treatment planning, immobilization, and plan-adaption criteria must be thoroughly explored to ensure adequate treatments are performed. This paper presents an institution's experience in developing a MRgART program for a 1.5T MR-linac for the first 234 patients. The paper suggests practical treatment workflows and considerations for treating with MRgART at different anatomical sites, including imaging guidelines, patient immobilization, adaptive workflows, and utilization of bolus.

Analysis of patient-specific quality assurance for Elekta Unity adaptive plans using statistical process control methodology.

The Elekta Unity MR-linac utilizes daily magnetic resonance imaging (MRI) for online plan adaptation. In the Unity workflow, adapt to position (ATP) and adapt to shape (ATS) treatment planning options are available which represent a virtual shift or full re-plan with contour adjustments respectively. Both techniques generate a new intensity modulated radiation therapy (IMRT) treatment plan while the patient lies on the treatment table and thus adapted plans cannot be measured prior to treatment delivery. A statistical process control methodology was used to analyze 512 patient-specific IMRT QA measurements performed on the MR-compatible SunNuclear ArcCheck with a gamma criterion of 3%/2 mm using global normalization and a 10% low dose threshold. The lower control limit (LCL) was determined from 68 IMRT reference plan measurements, and a one-sided process capability ratio ( C p , l ) was used to assess the pass rates from 432 measured ATP and 80 measured ATS plans. Further analysis was performed to assess differences between SBRT or conventional fractionation pass rates and to determine whether there was any correlation between the pass rates and plan complexity. The LCL of the reference plans was determined to be a gamma pass rate of 0.958, and the C p , l of the measured ATP plans and measured ATS plans were determined to be 1.403 and 0.940 for ATP and ATS plans, respectively, while a C p , l of 0.902 and 1.383 was found for SBRT and conventional fractionations respectively. For plan complexity, no correlation was found between modulation degree and gamma pass rate, but a statistically significant correlation was observed between the beam-averaged aperture area and gamma pass rate. All adaptive plans passed the TG-218 guidelines, but the ATS and SBRT plans tended to have a smaller beam-averaged aperture area with slightly lower gamma pass rates.

Diagnosing atmospheric communication of a sealed monitor chamber.

Daily output variations of up to ±2% were observed for a protracted time on a Varian TrueBeam® STx; these output variations were hypothesized to be the result of atmospheric communication of the sealed monitor chamber. Daily changes in output relative to baseline, measured with an ionization chamber array (DQA3) and the amorphous silicon flat panel detector (IDU) on the TrueBeam®, were compared with daily temperature-pressure corrections (PTP ) determined from sensors within the DQA3. Output measurements were performed using a Farmer® ionization chamber over a 5-hour period, during which there was controlled variation in the monitor chamber temperature. The root mean square difference between percentage output change from baseline measured with the DQA3 and IDU was 0.50% over all measurements. Over a 7-month retrospective review of daily changes in output and PTP , weak correlation (R2  = 0.30) was observed between output and PTP for the first 5 months; for the final 2 months, daily output changes were linearly correlated with changes in PTP , with a slope of 0.84 (R2  = 0.89). Ionization measurements corrected for ambient temperature and pressure during controlled heating and cooling of the monitor chamber differed from expected values for a sealed monitor chamber by up to 4.6%, but were consistent with expectation for an air-communicating monitor chamber within uncertainty (1.3%, k = 2). Following replacement of the depressurized monitor chamber, there has been no correlation between daily percentage change in output and PTP (R2  = 0.09). The utility of control charts is demonstrated for earlier identification of changes in the sensitivity of a sealed monitor chamber.

Commissioning of a 1.5T Elekta Unity MR-linac: A single institution experience.

MR image-guided radiotherapy has the potential to improve patient care, but integration of an MRI scanner with a linear accelerator adds complexity to the commissioning process. This work describes a single institution experience of commissioning an Elekta Unity MR-linac, including mechanical testing, MRI scanner commissioning, and dosimetric validation. Mechanical testing included multileaf collimator (MLC) positional accuracy, measurement of radiation isocenter diameter, and MR-to-MV coincidence. Key MRI tests included magnetic field homogeneity, geometric accuracy, image quality, and the accuracy of navigator-triggered imaging for motion management. Dosimetric validation consisted of comparison between measured and calculated PDDs and profiles, IMRT measurements, and end-to-end testing. Multileaf collimator positional accuracy was within 1.0 mm, the measured radiation isocenter walkout was 0.20 mm, and the coincidence between MR and MV isocenter was 1.06 mm, which is accounted for in the treatment planning system (TPS). For a 350-mm-diameter spherical volume, the peak-to-peak deviation of the magnetic field homogeneity was 4.44 ppm and the geometric distortion was 0.8 mm. All image quality metrics were within ACR recommendations. Navigator-triggered images showed a maximum deviation of 0.42, 0.75, and 3.0 mm in the target centroid location compared to the stationary target for a 20 mm motion at 10, 15, and 20 breaths per minute, respectively. TPS-calculated PDDs and profiles showed excellent agreement with measurement. The gamma passing rate for IMRT plans was 98.4 ± 1.1% (3%/ 2 mm) and end-to-end testing of adapted plans showed agreement within 0.4% between ion-chamber measurement and TPS calculation. All credentialing criteria were satisfied in an independent end-to-end test using an IROC MRgRT phantom.

MOSFET dosimeter characterization in MR-guided radiation therapy (MRgRT) Linac.

PURPOSE: With the increasing use of MR-guided radiation therapy (MRgRT), it becomes important to understand and explore accuracy of medical dosimeters in the presence of magnetic field. The purpose of this work is to characterize metal-oxide-semiconductor field-effect transistors (MOSFETs) in MRgRT systems at 0.345 T magnetic field strength. METHODS: A MOSFET dosimetry system, developed by Best Medical Canada for in-vivo patient dosimetry, was used to study various commissioning tests performed on a MRgRT system, MRIdian® Linac. We characterized the MOSFET dosimeter with different cable lengths by determining its calibration factor, monitor unit linearity, angular dependence, field size dependence, percentage depth dose (PDD) variation, output factor change, and intensity modulated radiation therapy quality assurance (IMRT QA) verification for several plans. MOSFET results were analyzed and compared with commissioning data and Monte Carlo calculations. RESULTS: MOSFET measurements were not found to be affected by the presence of 0.345 T magnetic field. Calibration factors were similar for different cable length dosimeters either placed at the parallel or perpendicular direction to the magnetic field, with variations of less than 2%. The detector showed good linearity (R2  = 0.999) for 100-600 MUs range. Output factor measurements were consistent with ionization chamber data within 2.2%. MOSFET PDD measurements were found to be within 1% for 1-15 cm depth range in comparison to ionization chamber. MOSFET normalized angular response matched thermoluminescent detector (TLD) response within 5.5%. The IMRT QA verification data for the MRgRT linac showed that the percentage difference between ionization chamber and MOSFET was 0.91%, 2.05%, and 2.63%, respectively for liver, spine, and mediastinum. CONCLUSION: MOSFET dosimeters are not affected by the 0.345 T magnetic field in MRgRT system. They showed physics parameters and performance comparable to TLD and ionization chamber; thus, they constitute an alternative to TLD for real-time in-vivo dosimetry in MRgRT procedures.

Characterization and longitudinal assessment of daily quality assurance for an MR-guided radiotherapy (MRgRT) linac.

PURPOSE: To describe and characterize daily machine quality assurance (QA) for an MR-guided radiotherapy (MRgRT) linac system, in addition to reporting a longitudinal assessment of the dosimetric and mechanical stability over a 7-month period of clinical operation. METHODS: Quality assurance procedures were developed to evaluate MR imaging/radiation isocenter, imaging and patient handling system, and linear accelerator stability. A longitudinal assessment was characterized for safety interlocks, laser and imaging isocenter coincidence, imaging and radiation (RT) isocentricity, radiation dose rate and output, couch motion, and MLC positioning. A cylindrical water phantom and an MR-compatible A1SL detector were utilized. MR and RT isocentricity and MLC positional accuracy was quantified through dose measured with a 0.40 cm2  x 0.83 cm2 field at each cardinal angle. The relationship between detector response to MR/RT isocentricity and MLC positioning was established through introducing known errors in phantom position. RESULTS: Correlation was found between detector response and introduced positional error (N = 27) with coefficients of determination of 0.9996 (IEC-X), 0.9967 (IEC-Y), 0.9968 (IEC-Z) in each respective shift direction. The relationship between dose (DoseMR/RT+MLC ) and the vector magnitude of MLC and MR/RT positional error (Errormag ) was calculated to be a nonlinear response and resembled a quadratic function: DoseMR/RT+MLC [%] = -0.0253 Errormag [mm]2  - 0.0195 Errormag [mm]. For the temporal assessment (N = 7 months), safety interlocks were functional. Laser coincidence to MR was within ±2.0 mm (99.6%) and ±1.0 mm (86.8%) over the 7-month assessment. IGRT position-reposition shifts were within ±2.0 mm (99.4%) and ±1.0 mm (92.4%). Output was within ±3% (99.4%). Mean MLC and MR/RT isocenter accuracy was 1.6 mm, averaged across cardinal angles for the 7-month period. CONCLUSIONS: The linac and IGRT accuracy of an MR-guided radiotherapy system has been validated and monitored over seven months for daily QA. Longitudinal assessment demonstrated a drift in dose rate, but temporal assessment of output, MLC position, and isocentricity has been stable.

Radiation treatment planning and delivery strategies for a pregnant brain tumor patient.

The management of a pregnant patient in radiation oncology is an infrequent event requiring careful consideration by both the physician and physicist. The aim of this manuscript was to highlight treatment planning techniques and detail measurements of fetal dose for a pregnant patient recently requiring treatment for a brain cancer. A 27-year-old woman was treated during gestational weeks 19-25 for a resected grade 3 astrocytoma to 50.4 Gy in 28 fractions, followed by an additional 9 Gy boost in five fractions. Four potential plans were developed for the patient: a 6 MV 3D-conformal treatment plan with enhanced dynamic wedges, a 6 MV step-and-shoot (SnS) intensity-modulated radiation therapy (IMRT) plan, an unflattened 6 MV SnS IMRT plan, and an Accuray TomoTherapy HDA helical IMRT treatment plan. All treatment plans used strategies to reduce peripheral dose. Fetal dose was estimated for each treatment plan using available literature references, and measurements were made using thermoluminescent dosimeters (TLDs) and an ionization chamber with an anthropomorphic phantom. TLD measurements from a full-course radiation delivery ranged from 1.0 to 1.6 cGy for the 3D-conformal treatment plan, from 1.0 to 1.5 cGy for the 6 MV SnS IMRT plan, from 0.6 to 1.0 cGy for the unflattened 6 MV SnS IMRT plan, and from 1.9 to 2.6 cGy for the TomoTherapy treatment plan. The unflattened 6 MV SnS IMRT treatment plan was selected for treatment for this particular patient, though the fetal doses from all treatment plans were deemed acceptable. The cumulative dose to the patient's unshielded fetus is estimated to be 1.0 cGy at most. The planning technique and distance between the treatment target and fetus both contributed to this relatively low fetal dose. Relevant treatment planning strategies and treatment delivery considerations are discussed to aid radiation oncologists and medical physicists in the management of pregnant patients.

A dynamic model-based approach to motion and deformation tracking of prosthetic valves from biplane x-ray images.

PURPOSE: Transcatheter aortic valve replacement (TAVR) is a minimally invasive procedure in which a prosthetic heart valve is placed and expanded within a defective aortic valve. The device placement is commonly performed using two-dimensional (2D) fluoroscopic imaging. Within this work, we propose a novel technique to track the motion and deformation of the prosthetic valve in three dimensions based on biplane fluoroscopic image sequences. METHODS: The tracking approach uses a parameterized point cloud model of the valve stent which can undergo rigid three-dimensional (3D) transformation and different modes of expansion. Rigid elements of the model are individually rotated and translated in three dimensions to approximate the motions of the stent. Tracking is performed using an iterative 2D-3D registration procedure which estimates the model parameters by minimizing the mean-squared image values at the positions of the forward-projected model points. Additionally, an initialization technique is proposed, which locates clusters of salient features to determine the initial position and orientation of the model. RESULTS: The proposed algorithms were evaluated based on simulations using a digital 4D CT phantom as well as experimentally acquired images of a prosthetic valve inside a chest phantom with anatomical background features. The target registration error was 0.12 ± 0.04 mm in the simulations and 0.64 ± 0.09 mm in the experimental data. CONCLUSIONS: The proposed algorithm could be used to generate 3D visualization of the prosthetic valve from two projections. In combination with soft-tissue sensitive-imaging techniques like transesophageal echocardiography, this technique could enable 3D image guidance during TAVR procedures.

Automated 3D coronary sinus catheter detection using a scanning-beam digital x-ray system.

Scanning-beam digital x-ray (SBDX) is an inverse geometry x-ray fluoroscopy system capable of tomosynthesis-based 3D tracking of catheter electrodes concurrent with fluoroscopic display. To facilitate respiratory motion-compensated 3D catheter tracking, an automated coronary sinus (CS) catheter detection algorithm for SBDX was developed. The technique uses the 3D localization capability of SBDX and prior knowledge of the catheter shape. Candidate groups of points representing the CS catheter are obtained from a 3D shape-constrained search. A cost function is then minimized over the groups to select the most probable CS catheter candidate. The algorithm was implemented in MATLAB and tested offline using recorded image sequences of a chest phantom containing a CS catheter, ablation catheter, and fiducial clutter. Fiducial placement was varied to create challenging detection scenarios. Table panning and elevation was used to simulate motion. The CS catheter detection method had 98.1% true positive rate and 100% true negative rate in 2755 frames of imaging. Average processing time was 12.7 ms/frame on a PC with a 3.4 GHz CPU and 8 GB memory. Motion compensation based on 3D CS catheter tracking was demonstrated in a moving chest phantom with a fixed CS catheter and an ablation catheter pulled along a fixed trajectory. The RMS error in the tracked ablation catheter trajectory was 1.41 mm, versus 10.35 mm without motion compensation. A computationally efficient method of automated 3D CS catheter detection has been developed to assist with motion-compensated 3D catheter tracking and registration of 3D cardiac models to tracked catheters.

Localization of cardiac volume and patient features in inverse geometry x-ray fluoroscopy.

The scanning-beam digital x-ray (SBDX) system is an inverse geometry x-ray fluoroscopy technology that performs real-time tomosynthesis at planes perpendicular to the source-detector axis. The live display is a composite image which portrays sharp features (e.g. coronary arteries) extracted from a 16 cm thick reconstruction volume. We present a method for automatically determining the position of the cardiac volume prior to acquisition of a coronary angiogram. In the algorithm, a single non-contrast frame is reconstructed over a 44 cm thickness using shift-and-add digital tomosynthesis. Gradient filtering is applied to each plane to emphasize features such as the cardiomediastinal contour, diaphragm, and lung texture, and then sharpness vs. plane position curves are generated. Three sharpness metrics were investigated: average gradient in the bright field, maximum gradient, and the number of normalized gradients exceeding 0.5. A model correlating the peak sharpness in a non-contrast frame and the midplane of the coronary arteries in a contrast-enhanced frame was established using 37 SBDX angiographic loops (64-136 kg human subjects, 0-30° cranial-caudal). The average gradient in the bright field (primarily lung) and the number of normalized gradients >0.5 each yielded peaks correlated to the coronary midplane. The rms deviation between the predicted and true midplane was 1.57 cm. For a 16 cm reconstruction volume and the 5.5-11.5 cm thick cardiac volumes in this study, midplane estimation errors of 2.25-5.25 cm were tolerable. Tomosynthesis-based localization of cardiac volume is feasible. This technique could be applied prior to coronary angiography, or to assist in isocentering the patient for rotational angiography.

Single-view geometric calibration for C-arm inverse geometry CT.

Accurate and artifact-free reconstruction of tomographic images requires precise knowledge of the imaging system geometry. A projection matrix-based calibration method to enable C-arm inverse geometry CT (IGCT) is proposed. The method is evaluated for scanning-beam digital x-ray (SBDX), a C-arm mounted inverse geometry fluoroscopic technology. A helical configuration of fiducials is imaged at each gantry angle in a rotational acquisition. For each gantry angle, digital tomosynthesis is performed at multiple planes and a composite image analogous to a cone-beam projection is generated from the plane stack. The geometry of the C-arm, source array, and detector array is determined at each angle by constructing a parameterized three-dimensional-to-two-dimensional projection matrix that minimizes the sum-of-squared deviations between measured and projected fiducial coordinates. Simulations were used to evaluate calibration performance with translations and rotations of the source and detector. The relative root-mean-square error in a reconstruction of a numerical thorax phantom was 0.4% using the calibration method versus 7.7% without calibration. In phantom studies, reconstruction of SBDX projections using the proposed method eliminated artifacts present in noncalibrated reconstructions. The proposed IGCT calibration method reduces image artifacts when uncertainties exist in system geometry.

Dynamic electronic collimation method for 3-D catheter tracking on a scanning-beam digital x-ray system.

Scanning-beam digital x-ray (SBDX) is an inverse geometry x-ray fluoroscopy system capable of tomosynthesis-based 3-D catheter tracking. This work proposes a method of dose-reduced 3-D catheter tracking using dynamic electronic collimation (DEC) of the SBDX scanning x-ray tube. This is achieved through the selective deactivation of focal spot positions not needed for the catheter tracking task. The technique was retrospectively evaluated with SBDX detector data recorded during a phantom study. DEC imaging of a catheter tip at isocenter required 340 active focal spots per frame versus 4473 spots in full field-of-view (FOV) mode. The dose-area product (DAP) and peak skin dose (PSD) for DEC versus full FOV scanning were calculated using an SBDX Monte Carlo simulation code. The average DAP was reduced to 7.8% of the full FOV value, consistent with the relative number of active focal spots (7.6%). For image sequences with a moving catheter, PSD was 33.6% to 34.8% of the full FOV value. The root-mean-squared-deviation between DEC-based 3-D tracking coordinates and full FOV 3-D tracking coordinates was less than 0.1 mm. The 3-D distance between the tracked tip and the sheath centerline averaged 0.75 mm. DEC is a feasible method for dose reduction during SBDX 3-D catheter tracking.

A geometric calibration method for inverse geometry computed tomography using P-matrices.

Accurate and artifact free reconstruction of tomographic images requires precise knowledge of the imaging system geometry. This work proposes a novel projection matrix (P-matrix) based calibration method to enable C-arm inverse geometry CT (IGCT). The method is evaluated for scanning-beam digital x-ray (SBDX), a C-arm mounted inverse geometry fluoroscopic technology. A helical configuration of fiducials is imaged at each gantry angle in a rotational acquisition. For each gantry angle, digital tomosynthesis is performed at multiple planes and a composite image analogous to a cone-beam projection is generated from the plane stack. The geometry of the C-arm, source array, and detector array is determined at each angle by constructing a parameterized 3D-to-2D projection matrix that minimizes the sum-of-squared deviations between measured and projected fiducial coordinates. Simulations were used to evaluate calibration performance with translations and rotations of the source and detector. In a geometry with 1 mm translation of the central ray relative to the axis-of-rotation and 1 degree yaw of the detector and source arrays, the maximum error in the recovered translational parameters was 0.4 mm and maximum error in the rotation parameter was 0.02 degrees. The relative root-mean-square error in a reconstruction of a numerical thorax phantom was 0.4% using the calibration method, versus 7.7% without calibration. Changes in source-detector-distance were the most challenging to estimate. Reconstruction of experimental SBDX data using the proposed method eliminated double contour artifacts present in a non-calibrated reconstruction. The proposed IGCT geometric calibration method reduces image artifacts when uncertainties exist in system geometry.

Method for dose-reduced 3D catheter tracking on a scanning-beam digital x-ray system using dynamic electronic collimation.

Scanning-beam digital x-ray (SBDX) is an inverse geometry x-ray fluoroscopy system capable of tomosynthesis-based 3D catheter tracking. This work proposes a method of dose-reduced 3D tracking using dynamic electronic collimation (DEC) of the SBDX scanning x-ray tube. Positions in the 2D focal spot array are selectively activated to create a region-of-interest (ROI) x-ray field around the tracked catheter. The ROI position is updated for each frame based on a motion vector calculated from the two most recent 3D tracking results. The technique was evaluated with SBDX data acquired as a catheter tip inside a chest phantom was pulled along a 3D trajectory. DEC scans were retrospectively generated from the detector images stored for each focal spot position. DEC imaging of a catheter tip in a volume measuring 11.4 cm across at isocenter required 340 active focal spots per frame, versus 4473 spots in full-FOV mode. The dose-area-product (DAP) and peak skin dose (PSD) for DEC versus full field-of-view (FOV) scanning were calculated using an SBDX Monte Carlo simulation code. DAP was reduced to 7.4% to 8.4% of the full-FOV value, consistent with the relative number of active focal spots (7.6%). For image sequences with a moving catheter, PSD was 33.6% to 34.8% of the full-FOV value. The root-mean-squared-deviation between DEC-based 3D tracking coordinates and full-FOV 3D tracking coordinates was less than 0.1 mm. The 3D distance between the tracked tip and the sheath centerline averaged 0.75 mm. Dynamic electronic collimation can reduce dose with minimal change in tracking performance.

Depth-resolved registration of transesophageal echo to x-ray fluoroscopy using an inverse geometry fluoroscopy system.

PURPOSE: Image registration between standard x-ray fluoroscopy and transesophageal echocardiography (TEE) has recently been proposed. Scanning-beam digital x-ray (SBDX) is an inverse geometry fluoroscopy system designed for cardiac procedures. This study presents a method for 3D registration of SBDX and TEE images based on the tomosynthesis and 3D tracking capabilities of SBDX. METHODS: The registration algorithm utilizes the stack of tomosynthetic planes produced by the SBDX system to estimate the physical 3D coordinates of salient key-points on the TEE probe. The key-points are used to arrive at an initial estimate of the probe pose, which is then refined using a 2D/3D registration method adapted for inverse geometry fluoroscopy. A phantom study was conducted to evaluate probe pose estimation accuracy relative to the ground truth, as defined by a set of coregistered fiducial markers. This experiment was conducted with varying probe poses and levels of signal difference-to-noise ratio (SDNR). Additional phantom and in vivo studies were performed to evaluate the correspondence of catheter tip positions in TEE and x-ray images following registration of the two modalities. RESULTS: Target registration error (TRE) was used to characterize both pose estimation and registration accuracy. In the study of pose estimation accuracy, successful pose estimates (3D TRE < 5.0 mm) were obtained in 97% of cases when the SDNR was 5.9 or higher in seven out of eight poses. Under these conditions, 3D TRE was 2.32 ± 1.88 mm, and 2D (projection) TRE was 1.61 ± 1.36 mm. Probe localization error along the source-detector axis was 0.87 ± 1.31 mm. For the in vivo experiments, mean 3D TRE ranged from 2.6 to 4.6 mm and mean 2D TRE ranged from 1.1 to 1.6 mm. Anatomy extracted from the echo images appeared well aligned when projected onto the SBDX images. CONCLUSIONS: Full 6 DOF image registration between SBDX and TEE is feasible and accurate to within 5 mm. Future studies will focus on real-time implementation and application-specific analysis.

Detector, collimator and real-time reconstructor for a new scanning-beam digital x-ray (SBDX) prototype.

Scanning-beam digital x-ray (SBDX) is an inverse geometry fluoroscopy system for low dose cardiac imaging. The use of a narrow scanned x-ray beam in SBDX reduces detected x-ray scatter and improves dose efficiency, however the tight beam collimation also limits the maximum achievable x-ray fluence. To increase the fluence available for imaging, we have constructed a new SBDX prototype with a wider x-ray beam, larger-area detector, and new real-time image reconstructor. Imaging is performed with a scanning source that generates 40,328 narrow overlapping projections from 71 × 71 focal spot positions for every 1/15 s scan period. A high speed 2-mm thick CdTe photon counting detector was constructed with 320×160 elements and 10.6 cm × 5.3 cm area (full readout every 1.28 µs), providing an 86% increase in area over the previous SBDX prototype. A matching multihole collimator was fabricated from layers of tungsten, brass, and lead, and a multi-GPU reconstructor was assembled to reconstruct the stream of captured detector images into full field-of-view images in real time. Thirty-two tomosynthetic planes spaced by 5 mm plus a multiplane composite image are produced for each scan frame. Noise equivalent quanta on the new SBDX prototype measured 63%-71% higher than the previous prototype. X-ray scatter fraction was 3.9-7.8% when imaging 23.3-32.6 cm acrylic phantoms, versus 2.3-4.2% with the previous prototype. Coronary angiographic imaging at 15 frame/s was successfully performed on the new SBDX prototype, with live display of either a multiplane composite or single plane image.

Feasibility of CT-based 3D anatomic mapping with a scanning-beam digital x-ray (SBDX) system.

This study investigates the feasibility of obtaining CT-derived 3D surfaces from data provided by the scanning-beam digital x-ray (SBDX) system. Simulated SBDX short-scan acquisitions of a Shepp-Logan and a thorax phantom containing a high contrast spherical volume were generated. 3D reconstructions were performed using a penalized weighted least squares method with total variation regularization (PWLS-TV), as well as a more efficient variant employing gridding of projection data to parallel rays (gPWLS-TV). Voxel noise, edge blurring, and surface accuracy were compared to gridded filtered back projection (gFBP). PWLS reconstruction of a noise-free reduced-size Shepp-Logan phantom had 1.4% rRMSE. In noisy gPWLS-TV reconstructions of a reduced-size thorax phantom, 99% of points on the segmented sphere perimeter were within 0.33, 0.47, and 0.70 mm of the ground truth, respectively, for fluences comparable to imaging through 18.0, 27.2, and 34.6 cm acrylic. Surface accuracies of gFBP and gPWLS-TV were similar at high fluences, while gPWLS-TV offered improvement at the lowest fluence. The gPWLS-TV voxel noise was reduced by 60% relative to gFBP, on average. High-contrast linespread functions measured 1.25 mm and 0.96 mm (FWHM) for gPWLS-TV and gFBP. In a simulation of gated and truncated projection data from a full-sized thorax, gPWLS-TV reconstruction yielded segmented surface points which were within 1.41 mm of ground truth. Results support the feasibility of 3D surface segmentation with SBDX. Further investigation of artifacts caused by data truncation and patient motion is warranted.

Monte Carlo simulation of inverse geometry x-ray fluoroscopy using a modified MC-GPU framework.

Scanning-Beam Digital X-ray (SBDX) is a technology for low-dose fluoroscopy that employs inverse geometry x-ray beam scanning. To assist with rapid modeling of inverse geometry x-ray systems, we have developed a Monte Carlo (MC) simulation tool based on the MC-GPU framework. MC-GPU version 1.3 was modified to implement a 2D array of focal spot positions on a plane, with individually adjustable x-ray outputs, each producing a narrow x-ray beam directed toward a stationary photon-counting detector array. Geometric accuracy and blurring behavior in tomosynthesis reconstructions were evaluated from simulated images of a 3D arrangement of spheres. The artifact spread function from simulation agreed with experiment to within 1.6% (rRMSD). Detected x-ray scatter fraction was simulated for two SBDX detector geometries and compared to experiments. For the current SBDX prototype (10.6 cm wide by 5.3 cm tall detector), x-ray scatter fraction measured 2.8-6.4% (18.6-31.5 cm acrylic, 100 kV), versus 2.1-4.5% in MC simulation. Experimental trends in scatter versus detector size and phantom thickness were observed in simulation. For dose evaluation, an anthropomorphic phantom was imaged using regular and regional adaptive exposure (RAE) scanning. The reduction in kerma-area-product resulting from RAE scanning was 45% in radiochromic film measurements, versus 46% in simulation. The integral kerma calculated from TLD measurement points within the phantom was 57% lower when using RAE, versus 61% lower in simulation. This MC tool may be used to estimate tomographic blur, detected scatter, and dose distributions when developing inverse geometry x-ray systems.